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Observation of a Strongly-Interacting Degenerate Fermi Gas of Atoms 

K. M. O'Hara, S. L. Hemmer, M. E. Gehm, S. R. Granade, and J. E. Thomas 
Physics Department, Duke University, Durham, North Carolina 27708-0305 

We report on the observation of a highly-degenerate, strongly-interacting Fermi gas of atoms. 
Fermionic 6 Li atoms in an optical trap are evaporatively cooled to degeneracy using a magnetic 
field to induce strong, resonant interactions. Upon abruptly releasing the cloud from the trap, 
the gas is observed to expand rapidly in the transverse direction while remaining nearly stationary 
in the axial. We interpret the expansion dynamics in terms of collisionless superfluid and colli- 
sional hydrodynamics. For the data taken at the longest evaporation times, we find that collisional 
hydrodynamics does not provide a satisfactory explanation, while superfluidity is plausible. 



As the fundamental constituents of matter are in- 
teracting fermions, the experimental study of strongly- 
interacting, degenerate Fermi gases will impact theories 
in fields from particle physics to materials science. Al- 
though the interactions between fermions are understood 
when they are weak (e.g., quantum electrodynamics), 
the treatment of very strong interactions requires the de- 
velopment of new theoretical approaches. To test these 
new approaches, there is a need for experimental systems 
with widely tunable interaction strengths, densities, and 
temperatures. Ultracold atomic Fermi gases have ex- 
actly these properties, and thus enable tests of calcula- 
tional techniques for fundamental systems ranging from 
quarks in nuclear matter to electrons in high tempera- 
ture superconductors g g |, g g g, g, g |, 0, 0. For 
this reason, a number of groups are developing meth- 
ods for creating and exploring ultracold atomic Fermi 

gases II P, 0, HJ 0- We re P° rt on the stud y 
of a strongly-interacting, degenerate Fermi gas. In con- 
trast to the isotropic expansion previously observed for 
a noninteracting degenerate Fermi gas |l^| , we observe 
anisotropic expansion when the gas is released from an 
optical trap. 

An exciting feature of strongly-interacting atomic 
Fermi gases is the possibility of high-temperature su- 
perfluids that are analogs of very high temperature su- 
perconductors H g, |l^, pT[ . Our experiments produce 
the conditions predicted for this type of superfluid tran- 
sition. Further, the anisotropic expansion we observe 
has been suggested as a signature of the onset of super- 
fluidity in a Fermi gas [Q. We interpret the observed 
anisotropic expansion in terms of both collisionless su- 
perfluid hydrodynamics Jl8| and a new form of collisional 
hydrodynamics. 



Strong, magnetically-tunable interactions are achieved 
in our experiments by employing a Fermi gas compris- 
ing a 50-50 mixture of the two lowest hyperfine states 
of 6 Li, i.e., the \F = 1/2, M = ±1/2) states in the 
low-magnetic-field basis. This mixture has a predicted 
broad Feshbach resonance near an applied magnetic field 
of 860 G ||, |C}, where the energy of a bound 6 Li- 6 Li 
molecular state is tuned into coincidence with the total 
energy of the colliding atoms. This enables the interac- 

~ [HJ |§. It has 
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tion strength to be widely varied [H 
also been suggested that interactions between fermions 
can be modified by immersion in a Bose gas p3[ |. Our 
experiments are performed at 910 G, where the zero en- 



ergy scattering length as is estimated to be 
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(a = 0.53 x 10 -8 cm) and the gas has strongly attractive 
interactions. Resonance superfluidity has been predicted 
to occur at this magnetic field for sufficiently low tem- 
peratures fLlf . 

In our experiments, 6 Li atoms are loaded from a 
magneto-optical trap into an ultrastable CO2 laser 
trap Jig] . The trap oscillation frequencies are uj z = 
2tt x (230 ± 20 Hz) for the axial (z) direction and 
uj± = 2tt x (6625 ± 50 Hz) for the transverse direc- 
tions. Rate equation pumping is used to produce the 
50-50 spin-mixture: A broad-band radio frequency (rf) 
field centered at 7.4 MHz is applied at a magnetic field 
of ~ 8 G, nulling the population difference according 
to An(t) — An(0) exp(— 2Rt), where R is the pump- 
ing rate. In our experiments, 2R = 600 sec -1 ; applying 
the rf field for t = 0.1 sec produces precise population 
balance. 

We achieve very low temperatures via rapid forced 
evaporation in the CO2 laser trap. In contrast to ex- 
periments which employ magnetic traps to achieve de- 
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generacy [jl^, [l3[ [l4|, 0, this approach has several 
natural advantages. First, we are able to evaporate both 
spin states to degeneracy at the desired magnetic field of 
910 G. As a result, the sample is never exposed to fields 
near 650 G where loss and heating are observed jl9], . 
Second, the evaporation process is identical for both spin 
states, thereby maintaining the initial spin balance as 
well as Fermi surface matching. Third, at this field, the 
collision cross section is extremely large and unitarity- 
limited, so that runaway evaporation is expected psf . 

Forced evaporation is achieved by lowering the power 
of the trapping laser, while maintaining the beam pro- 
file and angular alignment. The trap depth U is re- 
duced for 3.5 s according to the trajectory U{t) = 
U (1 + t/r)- 1A5 + U B H, where U B is a small offset. 
The value of t is taken to be 0.1 s, large compared to 
the time constant estimated for achieving degeneracy at 
910 G. With this choice, very high evaporation efficiency 
is achieved, yielding extremely low temperatures. 

After evaporation, the trap is adiabatically recom- 
pressed to full depth over 0.5 s and then held for 0.5 
s to ensure thermal equilibrium. While maintaining the 
applied magnetic field of 910 G, the gas is released from 
the trap and imaged at various times to observe the 
anisotropy. The CO2 laser power is extinguished in less 
than 1 /is with a rejection ratio of 2 x 10~ 5 [ ^5| . 

A CCD camera images the gas from a direction per- 
pendicular to the axial direction (z) of the trap and par- 
allel to the applied magnetic field direction (y). The 
small repulsive potential induced by the high-field mag- 
net is along the camera observation axis and does not 
affect the images. The remaining attractive potential 
has cylindrical symmetry and corresponds to a harmonic 
potential with an oscillation frequency for 6 Li of 20 Hz 
at 910 G. Resonant absorption imaging is performed on 
a cycling transition at a fixed high magnetic field by 
using a weak (J// sa t — 0.05), 20 /is probe laser pulse 
that is cr_ polarized with respect to the y-axis. Any 
residual a + component is rejected by an analyzer. At 
910 G, the transitions originating from the two occupied 
spin states are split by 70 MHz and are well resolved 
compared to the 3 MHz half-linewidth, permitting pre- 
cise determination of the column density and hence the 
number of atoms per state. The magnification is found 
to be 4.9 ± 0.15 by moving the axial position of the trap 
through 0.5 mm with a micrometer. The net system- 
atic error in the number measurement is estimated to 
be |^6) ■ The spatial resolution is estimated to be 

~ 4 (im by quadratically combining the effective pixel 
size, 13.0 /zm/4.9, with the aperture limited spatial res- 
olution of ~ 3 fim. 
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FIG. 1: False-color absorption images of a strongly interact- 
ing, degenerate Fermi gas as a function of time t after release 
from full trap depth for t — 0.1 — 2.0 ms, top to bottom. 
The axial width of the gas remains nearly stationary as the 
transverse width expands rapidly. 



Figure 1 shows images of the anisotropic expansion of 
the degenerate gas at various times t after release from 
full trap depth. The gas rapidly expands in the trans- 
verse direction (Fig. 2A) while remaining nearly station- 
ary in the axial direction (Fig. 2B) over a time period 
of 2.0 ms. In contrast to ballistic expansion, where the 
column density is cx 1/i 2 , the column density decreases 
only as 1/t for anisotropic expansion. Consequently, the 
signals are quite large even for long expansion times. 

One possible explanation of the observed anisotropy 
is provided by a recent theory of collisionless superfluid 
hydrodynamics |[8| . After release from the trap, the gas 
expands hydrodynamically due to the force from an effec- 
tive potential U c fi = ep + Umf, where £f(x) is the local 
Fermi energy and C/mf(x) is the mean field contribution. 
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FIG. 2: One-dimensional spatial distributions in the trans- 
verse (A) and axial (B) directions (red, 0.4 ms; blue, 1.0 
ms; green, 2.0 ms). The transverse distributions are shown 
fit with zero-temperature Thomas-Fermi distributions, while 
the axial are shown fit with gaussian distributions. 



In general, Umf cx a c ff n, where n is the spatial density 
and a c g is an effective scattering length. However, this 
theory is not rigorously applicable to our experiment, as 
it was derived for the dilute limit assuming a momentum- 
independent scattering length a e g = as- This assump- 
tion is only valid when fcp|as| < 1, where kp is the 
Fermi wavevector. By contrast, our experiments are 
performed in the intermediate density regime , where 
fc^lasl S> 1, and the interactions are unitarity-limited. 
We have therefore attempted to extend the theory in 
the context of a simple model. We make the assumption 
that unitarity limits a c g to ~ 1/kp. As n oc kp and 
cf( x ) = h 2 kp/(2M), we obtain J7 MF = f3ep(x), where 
j3 is a constant. This simple assumption is further justi- 
fied by more detailed calculations which show that /3 
is an important universal many-body parameter. With 



this assumption, U e g (x) = (1 + 0) ep{x). As oc n 2 / 3 , 
it then follows that U c s oc n 2 / 3 . 

For release from a harmonic trap and U c g oc n 1 , the 
hydrodynamic equations admit an exact solution JlJ, 



n(x, t) 



n (x/b x ,y/b y ,z/b z ) 



b x b v b z 



(1) 



Here, no(x) is the initial spatial distribution in the trap, 
and bi(t) are time dependent scaling parameters, which 
satisfy simple coupled differential equations with the ini- 
tial conditions, frj(0) = 1, 6j(0) = 0, where i — x,y, z. As 
the shape of the initial distribution (even with the mean 
field included) is determined by the trap potential, no is 
a function only of r' where uj 2 r' 2 = ui\{x 2 + y 2 ) + lo 2 z 2 



and u> 



a; i uj. 



a/3 



Hence, the initial radii of the den- 



sity distribution no are in the proportion a x (0) / a z (0) = 
A = uj z /ujj_. For our trap, A = 0.035. Then, during 
hydrodynamic expansion, the radii of the density distri- 
bution evolve according to 



c (i) = a x (0) b x (t) , a z (t) = a z (0) b z (t) . 



(2) 



We determine b x {t) — b y (t) and b z (t) from their evolu- 



For 7 = 2/3, b x 



b x 7/ % 



2/3 



tion equations 
and b z — uj 2 z bx^'^bz 

^From the expansion data, the widths a x {t) and 
cr z (t) are determined by fitting one-dimensional distri- 
butions p8| with normalized, zero-temperature Thomas- 
Fermi (T-F) distributions, n(x)/N = (16/(57rcr x ))(l - 
x 2 /a 2 ) 5 / 2 . As shown in Fig. 2A, the zero-temperature 
T-F fits to the transverse spatial profiles are quite good. 
This shape is not unreasonable despite a potentially large 
mean field interaction. As noted above !7mfW oc ep(pc). 
Hence, the mean field simply rescales the Fermi energy 
in the equation of state pq| . In this case, it is easy to 
show that the initial shape of the cloud is expected to 
be that of a Thomas-Fermi distribution. This shape is 
then maintained by the hydrodyamic scaling of Eq. |^. 

Figure 3A shows the measured values of a x (t) and 
<r z (t) as a function of time t after release. To com- 
pare these results with the predictions of Eq. ^, we take 
the initial dimensions of the cloud, cr^O) and cr z (0), 
to be the zero-temperature Fermi radii. For our mea- 
sured number N — 7.5j^o 5 x 10 4 atoms per state, and 
uj = 2tt x (2160 ± 65 Hz), the Fermi temperature is 
Tp = hb>{<c,N) l l*/kB = 7.9tg;!/iK at full trap depth. 
One then obtains cr x (0) = y / 2k B T F /Muj 2 = 3.6±0.1 ^m 
in the transverse direction, and cr z (0) = 103 ± 3/im in 
the axial. For these initial dimensions, we obtain very 



good agreement with our measurements using no free 
parameters, as shown by the solid curves in Fig. 3A. 

Figure 3B shows the measured aspect ratios 
<?x (t) I a z (t) and the theoretical predictions based on hy- 
drodynamic, ballistic, and attractive (f3 = —0.4) or re- 
pulsive (P — 0.4) collisionless mean field scaling [ |l8| . 
The observed expansion appears to be nearly hydro- 
dynamic. For comparison, we also show the mea- 
sured aspect ratios obtained for release at 530 G, where 
the scattering length has been measured to be nearly 
zero pl| . In this case, there is excellent agree- 
ment with the ballistic expansion expected for a nonin- 
teracting gas. This directly confirms that the observed 
anisotropy is a consequence of interactions. 

A primary assumption of our simple model is that the 
system is unitarity limited, i.e., \kp as\ S> 1. In this limit 
the fluid properties are predicted to be independent of 
kpas, for both a collisionless superfluid [Q and for col- 
lisional hydrodynamics as shown below (see Eq. [|). To 
test this assumption, we have released the gas at a trap 
depth U which is 1/100 of the full trap depth. As the 
trap frequencies, and hence the Fermi energy, scale as 
VU, k F is reduced by (1/100) 1 ' A ~ 1/3. The time scale 
for the expansion is increased by a factor of 10 as ex- 
pected, and we continue to observe strongly anisotropic 
expansion. This suggests that the system is unitarity- 
limitcd at full trap depth, consistent with the estimated 
value of fej?|as| = 7.4 for our experimental conditions. 
Hence, the gas appears to be strongly interacting. 

In the unitarity limit, where Umf = @£f, we can ob- 
tain the first estimate of (3 from measurements of the 
transverse release energy. From the zero-temperature 
T-F fits to the 910 G data, we obtain an average trans- 
verse release energy of e x — (Mt)^/2) = M(a x /4:t) 2 = 
1.44 ± 0.02 /iK per particle. To derive an expression for 
the release energy as a function of /3, we assume that 
the initial density distribution no(x) is determined by 
the equation of state for a trapped, normal fluid | fl8| : 
£f(x) + PmfW + f/ t rap(x) = [i, where \i is the chemical 
potential. In this case, the release energy per particle is 

e r = | k B T F ^TTp. (3) 

As shown above, Tp = 7.91q2 at full trap depth. To 
obtain an estimate of (3 from the experiments, we as- 
sume e r /2 is released in each transverse direction, with 
negligible energy deposited axially. Calculating (3 us- 
ing the appropriate Tp and e x for each trial in the range 
t = 0.4-0.8 ms, we find = -0.10±0.07. Note that the 
uncertainty is the quadratic combination of the statisti- 
cal and systematic uncertainties in the measurements, 
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FIG. 3: (A) Transverse (red) and axial (blue) widths as func- 
tions of time after release. The solid curves are theoretical 
predictions based on hydrodynamic scaling with no free pa- 
rameters. (B) Aspect ratio of the cloud as a function of time 
after release. The dots indicate experimental data and the 
solid curves show theoretical predictions with no adjustable 
parameters (red, hydrodynamic; blue, ballistic; green, attrac- 
tive mean field; orange, repulsive mean field). 

but does not reflect any systematic effects arising from 
our model. We find that the sign is in agreement with 
recent predictions, but the magnitude is a factor of ten 
smaller than expected 0. It is remarkable that the re- 
lease energy and the initial cloud dimensions are well- 
described assuming a zero-temperature, noninteracting 
Fermi gas, while the highly anisotropic expansion results 
from strong interactions. 

The preceding determination of [3 assumed a cloud at 
zero temperature. To investigate the validity of this ap- 
proximation, we fit normalized, finite-temperature T-F 
distributions to the transverse data at short times, 0.2 
and 0.4 ms, where the signal to background ratio in the 
thermal wings is high. As noted above, it is not unrea- 
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sonable to expect a T-F distribution. In this case, it 
is possible to perform a two parameter fit with a x and 
T/Tp as the free parameters, where T/Tp is the ratio of 
the temperature to the Fermi temperature. The fits (not 
shown) to the measured transverse spatial distributions 
yield 0.08 < T/T F < 0.18, while a x is essentially un- 
changed from the zero-temperature results. Hence, it ap- 
pears that T/Tp is quite small and the zero-temperature 
approximation is reasonable. 

In our experiments, the peak Fermi density is calcu- 
lated to be n F = 4.7 x 10 13 /cm 3 per state |||. At 
this density, it is possible that the anisotropic expan- 
sion arises from collisional hydrodynamics [^J. In con- 
trast to the case usually considered, for the large scat- 
tering lengths in this system, the collision cross section 
is unitarity-limited to a value ap — An/kp in the de- 
generate regime. The gas is collisionally hydrodynamic 
when the collision parameter (j) = 7 /lo±_ 1, where 7 is 
the elastic collision rate and 1/uj± is the relevant time 
scale for the expansion. Using the s-wave Boltzmann 
equation p0| for a constant cross section a and includ- 
ing Pauli blocking p], we find that 7 = j F P {T/TF), 
where 70 = NMoui^/I/It: 2 kpTp), N is number of atoms 
in one spin state, and Fp describes the temperature 
dependence. Here, Fp — > T F /T for T > 2>, and 
F P ~ 15 (T/Tp) 2 for T < 0.2 Tp, where Pauli block- 
ing occurs. The maximum value of Fp = 1.3 occurs at 
T/Tp = 0.5. In the degenerate regime (i.e., T < 0.5 Tp), 
where a ~ ap, we obtain 



(6AJV) 



1/3 



Gtt 



-F P {T/Tp). 



(4) 



For our experimental conditions, (f> ~ 1.3 F P (T/Tp), 
which is independent of the trap depth as long as 
&f|os| ^ 1- As Fp is at most of order unity, strongly 
hydrodynamic behavior arising from collisions seems un- 
likely. Including the temperature dependence, 4> ranges 
from 0.8 down to 0.2 where the system is nearly collision- 
less. Hence, collisional hydrodynamics does not provide 
a satisfactory explanation of the observed anisotropic ex- 
pansion, while superfluid hydrodynamics is plausible. 

Given this possibility, we have performed an initial 
investigation of the transition between ballistic and hy- 
drodynamic expansion. We measure the aspect ratio for 
an expansion time of 0.6 ms as a function of the evapo- 
ration time. For short evaporation times < 0.13 s, where 
T/Tp > 3.5, the measured aspect ratio is consistent with 
that expected for ballistic expansion. For any evapora- 
tion time > 1.5 s, the aspect ratio is consistent with 
hydrodynamic expansion. We observe a very smooth 
transition between these two extremes. In the interme- 



diate regime, at temperatures below T/Tp = 3.5, the 
expansion lies between hydrodynamic and ballistic. At 
T/Tp = 3.5, where the evaporation time is short and 
the number is large, an estimate of the classical col- 
lision rate with a unitarity-limited cross section shows 
that the onset of collisional behavior is not surprising. 
In the intermediate region, there is no theory of expan- 
sion to describe the spatial anisotropy of the energy re- 
lease. Hence, any attempt to determine the tempera- 
ture is highly model-dependent and cannot be trusted. 
To further complicate the analysis, varying the evapora- 
tion time changes the trap population in addition to the 
temperature. Finally, if high-temperature resonance su- 
perfluidity does exist, the transition temperature is pre- 
dicted to be in the range 0.25-0.5 Tp, where Pauli block- 
ing is not very effective. Hence, one would not expect 
to observe a collisionless region immediately prior to the 
onset of superfluid hydrodynamics, unless the transition 
occurs at very low temperature, in contrast to predic- 
tions. 

There are a number of noticeable discrepancies be- 
tween the hydrodynamic theory and the data. The de- 
viations at and 0.1 ms can be explained by possible 
index of refraction effects as well as spatial resolution 
limits. These issues are not significant for longer expan- 
sion times where the density is reduced and the cloud size 
is well beyond the resolution limit of our imaging system. 
However, close examination of the long time deviations 
reveals that there may be a two-component structure 
in the gas. In the axial direction, hydrodynamic ex- 
pansion is very slow, and a second component expand- 
ing according to ballistic or collisionless mean field scal- 
ing (Fig. 3B) easily overtakes the hydrodynamic compo- 
nent. A two-component structure may also explain why 
the axial spatial distributions (Fig. 2B), are better fit 
by gaussian distributions than by zero-temperature T-F 
distributions. By contrast, in the transverse direction, 
the hydrodynamic expansion is the fastest, masking any 
two-component structure after a short time. 
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